home *** CD-ROM | disk | FTP | other *** search
/ The Datafile PD-CD 1 Issue 2 / PDCD-1 - Issue 02.iso / _utilities / utilities / 003 / sox_2 / sox / c / au < prev    next >
Text File  |  1993-09-23  |  8KB  |  318 lines
  1. /*
  2.  * Copyright 1991, 1992, 1993 Guido van Rossum And Sundry Contributors.
  3.  * This source code is freely redistributable and may be used for
  4.  * any purpose.  This copyright notice must be maintained. 
  5.  * Guido van Rossum And Sundry Contributors are not responsible for 
  6.  * the consequences of using this software.
  7.  */
  8.  
  9. /*
  10.  * Sound Tools Sun format with header (SunOS 4.1; see /usr/demo/SOUND).
  11.  * NeXT uses this format also, but has more format codes defined.
  12.  * DEC uses a slight variation and swaps bytes.
  13.  * We only support the common formats.
  14.  * CCITT G.721 (32 kbit/s) and G.723 (24/40 kbit/s) are also supported,
  15.  * courtesy Sun's public domain implementation.
  16.  * Output is always in big-endian (Sun/NeXT) order.
  17.  */
  18.  
  19. #include "st.h"
  20. #include "g72x.h"
  21.  
  22. /* Magic numbers used in Sun and NeXT audio files */
  23. #define SUN_MAGIC     0x2e736e64        /* Really '.snd' */
  24. #define SUN_INV_MAGIC    0x646e732e        /* '.snd' upside-down */
  25. #define DEC_MAGIC    0x2e736400        /* Really '\0ds.' (for DEC) */
  26. #define DEC_INV_MAGIC    0x0064732e        /* '\0ds.' upside-down */
  27. #define SUN_HDRSIZE    24            /* Size of minimal header */
  28. #define SUN_UNSPEC    ((unsigned)(~0))    /* Unspecified data size */
  29. #define SUN_ULAW    1            /* u-law encoding */
  30. #define SUN_LIN_8    2            /* Linear 8 bits */
  31. #define SUN_LIN_16    3            /* Linear 16 bits */
  32. #define SUN_LIN_24    4            /* Linear 24 bits */
  33. #define SUN_LIN_32    5            /* Linear 32 bits */
  34. #define SUN_FLOAT    6            /* IEEE FP 32 bits */
  35. #define SUN_DOUBLE    7            /* IEEE FP 64 bits */
  36. #define SUN_G721    23            /* CCITT G.721 4-bits ADPCM */
  37. #define SUN_G723_3    25            /* CCITT G.723 3-bits ADPCM */
  38. #define SUN_G723_5    26            /* CCITT G.723 5-bits ADPCM */
  39. /* The other formats are not supported by sox at the moment */
  40.  
  41. /* Private data */
  42. struct aupriv {
  43.     /* For writer: */
  44.     unsigned long data_size;
  45.     /* For G72x decoding: */
  46.     struct g72x_state state;
  47.     int (*dec_routine)();
  48.     int dec_bits;
  49.     unsigned int in_buffer;
  50.     int in_bits;
  51. };
  52.  
  53. IMPORT auwriteheader(P2(ft_t ft, unsigned long data_size));
  54.  
  55. austartread(ft) 
  56. ft_t ft;
  57. {
  58.     /* The following 6 variables represent a Sun sound header on disk.
  59.        The numbers are written as big-endians.
  60.        Any extra bytes (totalling hdr_size - 24) are an
  61.        "info" field of unspecified nature, usually a string.
  62.        By convention the header size is a multiple of 4. */
  63.     unsigned long magic;
  64.     unsigned long hdr_size;
  65.     unsigned long data_size;
  66.     unsigned long encoding;
  67.     unsigned long sample_rate;
  68.     unsigned long channels;
  69.  
  70.     register int i;
  71.     char *buf;
  72.     struct aupriv *p = (struct aupriv *) ft->priv;
  73.  
  74.     /* Sanity check */
  75.     if (sizeof(struct aupriv) > PRIVSIZE)
  76.         fail(
  77. "struct aupriv is too big (%d); change PRIVSIZE in st.h and recompile sox",
  78.              sizeof(struct aupriv));
  79.  
  80.     /* Check the magic word */
  81.     magic = rlong(ft);
  82.     if (magic == DEC_INV_MAGIC) {
  83.         ft->swap = 1;
  84.         report("Found inverted DEC magic word");
  85.     }
  86.     else if (magic == SUN_INV_MAGIC) {
  87.         ft->swap = 1;
  88.         report("Found inverted Sun/NeXT magic word");
  89.     }
  90.     else if (magic == SUN_MAGIC) {
  91.         ft->swap = 0;
  92.         report("Found Sun/NeXT magic word");
  93.     }
  94.     else if (magic == DEC_MAGIC) {
  95.         ft->swap = 0;
  96.         report("Found DEC magic word");
  97.     }
  98.     else
  99.         fail("Sun/NeXT/DEC header doesn't start with magic word\nTry the '.ul' file type with '-t ul -r 8000 filename'");
  100.  
  101.     /* Read the header size */
  102.     hdr_size = rlong(ft);
  103.     if (hdr_size < SUN_HDRSIZE)
  104.         fail("Sun/NeXT header size too small.");
  105.  
  106.     /* Read the data size; may be ~0 meaning unspecified */
  107.     data_size = rlong(ft);
  108.  
  109.     /* Read the encoding; there are some more possibilities */
  110.     encoding = rlong(ft);
  111.  
  112.  
  113.     /* Translate the encoding into style and size parameters */
  114.     /* (Or, for G.72x, set the decoding routine and parameters) */
  115.     p->dec_routine = NULL;
  116.     p->in_buffer = 0;
  117.     p->in_bits = 0;
  118.     switch (encoding) {
  119.     case SUN_ULAW:
  120.         ft->info.style = ULAW;
  121.         ft->info.size = BYTE;
  122.         break;
  123.     case SUN_LIN_8:
  124.         ft->info.style = SIGN2;
  125.         ft->info.size = BYTE;
  126.         break;
  127.     case SUN_LIN_16:
  128.         ft->info.style = SIGN2;
  129.         ft->info.size = WORD;
  130.         break;
  131.     case SUN_G721:
  132.         ft->info.style = SIGN2;
  133.         ft->info.size = WORD;
  134.         g72x_init_state(&p->state);
  135.         p->dec_routine = g721_decoder;
  136.         p->dec_bits = 4;
  137.         break;
  138.     case SUN_G723_3:
  139.         ft->info.style = SIGN2;
  140.         ft->info.size = WORD;
  141.         g72x_init_state(&p->state);
  142.         p->dec_routine = g723_24_decoder;
  143.         p->dec_bits = 4;
  144.         break;
  145.     case SUN_G723_5:
  146.         ft->info.style = SIGN2;
  147.         ft->info.size = WORD;
  148.         g72x_init_state(&p->state);
  149.         p->dec_routine = g723_40_decoder;
  150.         p->dec_bits = 4;
  151.         break;
  152.     default:
  153.         report("encoding: 0x%lx", encoding);
  154.         fail("Unsupported encoding in Sun/NeXT header.\nOnly U-law, signed bytes, and signed words are supported.");
  155.         /*NOTREACHED*/
  156.     }
  157.  
  158.     /* Read the sampling rate */
  159.     sample_rate = rlong(ft);
  160.     ft->info.rate = sample_rate;
  161.  
  162.     /* Read the number of channels */
  163.     channels = rlong(ft);
  164.     ft->info.channels = channels;
  165.  
  166.     /* Skip the info string in header; print it if verbose */
  167.     hdr_size -= SUN_HDRSIZE; /* #bytes already read */
  168.     if (hdr_size > 0) {
  169.         buf = (char *) malloc(hdr_size + 1);
  170.         for(i = 0; i < hdr_size; i++) {
  171.             buf[i] = (char) getc(ft->fp);
  172.             if (feof(ft->fp))
  173.                 fail("Unexpected EOF in Sun/NeXT header info.");
  174.         }
  175.         buf[i] = '\0';
  176.         ft->comment = buf;
  177.         report("Input file %s: Sun header info: %s", ft->filename, buf);
  178.     }
  179. }
  180.  
  181. /* When writing, the header is supposed to contain the number of
  182.    data bytes written, unless it is written to a pipe.
  183.    Since we don't know how many bytes will follow until we're done,
  184.    we first write the header with an unspecified number of bytes,
  185.    and at the end we rewind the file and write the header again
  186.    with the right size.  This only works if the file is seekable;
  187.    if it is not, the unspecified size remains in the header
  188.    (this is legal). */
  189.  
  190. austartwrite(ft) 
  191. ft_t ft;
  192. {
  193.     struct aupriv *p = (struct aupriv *) ft->priv;
  194.     int littlendian = 0;
  195.     char *endptr;
  196.  
  197.     p->data_size = 0;
  198.     auwriteheader(ft, SUN_UNSPEC);
  199.     endptr = (char *) &littlendian;
  200.     *endptr = 1;
  201.     if (littlendian == 1)
  202.         ft->swap = 1;
  203. }
  204.  
  205. /*
  206.  * Unpack input codes and pass them back as bytes.
  207.  * Returns 1 if there is residual input, returns -1 if eof, else returns 0.
  208.  * (Adapted from Sun's decode.c.)
  209.  */
  210. int
  211. unpack_input(ft, code)
  212. ft_t            ft;
  213. unsigned char        *code;
  214. {
  215.     struct aupriv        *p = (struct aupriv *) ft->priv;
  216.     unsigned char        in_byte;
  217.  
  218.     if (p->in_bits < p->dec_bits) {
  219.         if (fread(&in_byte, sizeof (char), 1, ft->fp) != 1) {
  220.             *code = 0;
  221.             return (-1);
  222.         }
  223.         p->in_buffer |= (in_byte << p->in_bits);
  224.         p->in_bits += 8;
  225.     }
  226.     *code = p->in_buffer & ((1 << p->dec_bits) - 1);
  227.     p->in_buffer >>= p->dec_bits;
  228.     p->in_bits -= p->dec_bits;
  229.     return (p->in_bits > 0);
  230. }
  231.  
  232. auread(ft, buf, samp)
  233. ft_t ft;
  234. long *buf, samp;
  235. {
  236.     struct aupriv *p = (struct aupriv *) ft->priv;
  237.     unsigned char code;
  238.     int done;
  239.     if (p->dec_routine == NULL)
  240.         return rawread(ft, buf, samp);
  241.     done = 0;
  242.     while (samp > 0 && unpack_input(ft, &code) >= 0) {
  243.         *buf++ = LEFT((*p->dec_routine)(code, AUDIO_ENCODING_LINEAR,
  244.                         &p->state),
  245.                   16);
  246.         samp--;
  247.         done++;
  248.     }
  249.     return done;
  250. }
  251.  
  252. auwrite(ft, buf, samp)
  253. ft_t ft;
  254. long *buf, samp;
  255. {
  256.     struct aupriv *p = (struct aupriv *) ft->priv;
  257.     p->data_size += samp * ft->info.size;
  258.     rawwrite(ft, buf, samp);
  259. }
  260.  
  261. void
  262. austopwrite(ft)
  263. ft_t ft;
  264. {
  265.     struct aupriv *p = (struct aupriv *) ft->priv;
  266.     if (!ft->seekable)
  267.         return;
  268.     if (fseek(ft->fp, 0L, 0) != 0)
  269.         fail("Can't rewind output file to rewrite Sun header.");
  270.     auwriteheader(ft, p->data_size);
  271. }
  272.  
  273. auwriteheader(ft, data_size)
  274. ft_t ft;
  275. unsigned long data_size;
  276. {
  277.     unsigned long magic;
  278.     unsigned long hdr_size;
  279.     unsigned long encoding;
  280.     unsigned long sample_rate;
  281.     unsigned long channels;
  282.  
  283.     if (ft->info.style == ULAW && ft->info.size == BYTE)
  284.         encoding = SUN_ULAW;
  285.     else if (ft->info.style == SIGN2 && ft->info.size == BYTE)
  286.         encoding = SUN_LIN_8;
  287.     else if (ft->info.style == SIGN2 && ft->info.size == WORD)
  288.         encoding = SUN_LIN_16;
  289.     else {
  290. /*        fail("Unsupported output style/size for Sun/NeXT header.  \nOnly U-law, signed bytes, and signed words are supported."); /* */
  291.         report("SUN .AU format not specified.  Defaulting to 8khz u-law");
  292.         encoding = SUN_ULAW;
  293.         ft->info.style = ULAW;
  294.         ft->info.size = BYTE;
  295.         ft->info.rate = 8012;  /* strange but true */
  296.     }
  297.  
  298.     magic = SUN_MAGIC;
  299.     wblong(ft, magic);
  300.  
  301.     /* hdr_size = SUN_HDRSIZE;        /* + strlen(ft->comment); */
  302.     hdr_size = SUN_HDRSIZE + strlen(ft->comment);
  303.     wblong(ft, hdr_size);
  304.  
  305.     wblong(ft, data_size);
  306.  
  307.     wblong(ft, encoding);
  308.  
  309.     sample_rate = ft->info.rate;
  310.     wblong(ft, sample_rate);
  311.  
  312.     channels = ft->info.channels;
  313.     wblong(ft, channels);
  314.  
  315.     fputs(ft->comment, ft->fp);
  316. }
  317.  
  318.